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Abstract 

We study the dual cascade scenario for two-dimensional turbulence driven by a 
spectrally localized forcing applied over a finite wavenumber range [&min; &max] (with 
^min > 0) such that the respective energy and enstrophy injection rates e and 
r) satisfy A^ in e < rj < £^ ax e. The classical Kraichnan-Leith-Batchelor paradigm, 
based on the simultaneous conservation of energy and enstrophy and the scale- 
selectivity of the molecular viscosity, requires that the domain be unbounded in 
both directions. For two-dimensional turbulence either in a doubly periodic domain 
or in an unbounded channel with a periodic boundary condition in the across- 
channel direction, a direct enstrophy cascade is not possible. In the usual case where 
the forcing wavenumber is no greater than the geometric mean of the integral and 
dissipation wavenumbers, constant spectral slopes must satisfy (3 > 5 and a + 
(5 > 8, where —a {—(5) is the asymptotic slope of the range of wavenumbers lower 
(higher) than the forcing wavenumber. The influence of a large-scale dissipation 
on the realizability of a dual cascade is analyzed. We discuss the consequences for 
numerical simulations attempting to mimic the classical unbounded picture in a 
bounded domain. 
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1 Introduction 



Since Kraichnan [13,14], Batchelor [2], and Leith [15] (referred to as KLB) 
adapted Kolmogorov's theory of self-similarity in three-dimensional turbulence 
to two-dimensional (2D) fluids, the conventional wisdom for decades has been 
that 2D turbulence simultaneously exhibits a direct cascade of enstrophy to 
large wavenumbers, up to a dissipation wavenumber k u , and an inverse cascade 
of energy to small wavenumbers, down to wavenumber k = 0. In the limit 
of small viscosity, the inverse cascade is thought to proceed indefinitely in 
time to ever-larger scales, transferring virtually all of the energy input to 
wavenumber zero. The direct cascade is thought to come into balance with 
viscosity, transferring virtually all of the enstrophy input to k u , where it will 
then be dissipated. In the long-time limit, a quasi-steady state is reached, in 
which two inertial ranges are established. According to the KLB theory, the 
energy range, which is only quasi-steady, scales as k~ 5 ^ 3 . The enstrophy range, 
which is in absolute equilibrium, should scale as k~ 3 . 

The idea of a dual cascade was first suggested by Fj0rtoft [10], who examined 
nonlinear transfer by individual interacting wavenumber triads. A later study 
by Merilees and Warn [23] provides more quantitative detail. They showed that 
roughly 70% (60%) of triads containing a given intermediate wavenumber k 
predominantly exchange energy (enstrophy) with lower (higher) wavenumbers. 
These analyses of nonlinear transfer, although carried out at the triad level, 
provide only a necessary basis for the KLB theory. The theory requires all 
triads to work collectively such that virtually all of the energy (enstrophy) 
input gets transferred to the large (small) scales; this is neither suggested nor 
implied by [10,23]. Nevertheless, a dual cascade from a spectrally localized 
initial spectrum is consistent with [10,23] and has been well confirmed by 
numerical simulations with various resolutions (e.g. Borue [4,5]; Frisch and 
Sulem [12]; Lilly [17]; Smith and Yakhot [34]). In particular, the inverse energy 
cascade is observed in the laboratory experiments of Dubos et al. [8] and Paret 
and Tabeling [26,27] What has not been established beyond doubt is the 
realization of the inertial spectral scaling k~ 3 . In fact, there exist other theories 
that propose very different spectral slopes (Moffatt [24]; Saffman [31]; Sulem 
and Frisch [37]). 

Numerical simulations, aiming to verify the KLB picture, face a number of 
formidable tasks. First, the simulations are performed for fluids in a rectan- 
gular box instead of an infinite domain. This turns out to be a serious short- 



1 The inverse energy cascade and the /c -5 / 3 range are also seen to be robust in 
many numerical simulations, at least until the energy reaches the lower spectral 
boundary; however, the evidence for a direct enstrophy cascade is inconclusive (see 
the discussion in Paret and Tabeling [27] and references therein). 



2 



coming, as the equilibrium dynamics of a fluid in a doubly periodic domain 
differs considerably from that of an unbounded fluid in the KLB picture (see 
Constantin, Foias, and Manley [7]; Tran and Shepherd [39]). In particular, 
the k" 3 range and the direct enstrophy cascade have been shown to be unreal- 
izable in a bounded domain. The main reason seems to be that in the bounded 
case, where an absolute equilibrium will be reached, there is no analogue of 
a persistent upscale flow of energy that eventually evades viscous dissipation 
altogether. To mimic the inverse cascade, one needs to introduce a dissipation 
that removes energy at the large scales (cf. [36]). Linear Ekman drag (propor- 
tional to k° and restricted to a few low wavenumbers) and inverse viscosity 
(i> t j - k 2fl , with n < 0) have often been used for that purpose. Unfortunately, 
the inverse energy cascade, subject to this large-scale dissipation, carries with 
it a significant fraction of the enstrophy (in contrast to the asymptotic KLB 
inverse cascade, which carries no enstrophy). Moreover, there is no guarantee 
that the large-scale dissipation will absorb virtually all of the energy input. 
This is crucial in the KLB picture, as any fraction of the energy input that 
gets reflected would ultimately be trapped in the inertial ranges. The trapped 
energy, being in a virtually inviscid region, would then considerably change 
the dynamics and the slopes of the inertial ranges. Second, testing the theory 
requires the achievement of high Reynolds numbers, but current computers 
are only able to resolve relatively low Reynolds numbers. To overcome these 
resolution limitations, researchers often resort to introducing a hyperviscosity 
z/^/c 2/i , where the degree of viscosity // is greater than one. This numerical 
device helps to compress the dissipation range, allowing simulations to be 
performed at relatively low resolutions. It is hoped that the effect of this mod- 
ification on the inertial-range dynamics is negligible. However, a significant 
dissipation of energy at the small scales is inevitable in all numerical schemes. 

In this study, we establish a theoretical basis highlighting the dynamical dif- 
ferences between bounded and unbounded 2D turbulence. The consequences 
for numerical simulations that aim to verify the KLB theory will be addressed. 
In particular, we argue that instead of the familiar inertial-range spectral scal- 
ings k~ 5 ^ 3 and k~ 3 conjectured by KLB, one would expect steeper scalings, 
such as k~ 3 and k~ 5 , respectively, in these ranges for bounded systems in 
equilibrium. In fact, it is shown in [39] that the enstrophy range must be 
steeper than k~ 5 ; the physical implication is that no direct enstrophy cascade 
is possible. This result is easily generalized, using the Poincare inequality, 
to turbulence in semi- unbounded 2D fluids, i.e. fluids confined to unbounded 
channels with a periodic boundary condition in the across-channel direction. 
Moreover, we investigate the effectiveness of using a large-scale dissipation to 
obtain a dual cascade in a bounded fluid. This information should be partic- 
ularly useful for numerical simulations. 

In Section 2, we summarize the KLB theory and contrast it to the dynamics 
of a bounded fluid. We also show that the theory does not apply to a fluid in 
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an infinite channel with a periodic boundary condition in the across-channel 
direction. In Section 3, we review a result in [39] related to the unrealizability 
of a direct enstrophy cascade and discuss the fundamental differences between 
the dynamics of bounded fluids in equilibrium and that of unbounded fluids 
(or bounded fluids in transient phase) in the KLB picture. In Section 4, we 
derive a new constraint on the spectral slope of the energy range for bounded 
fluids in equilibrium and a condition for a persistent inverse energy cascade for 
unbounded fluids. In Section 5, we point out some effects of including a large- 
scale dissipation and the implications for numerical simulations attempting to 
verify the KLB picture. We conclude with some remarks in the final section. 
Further estimates on spectral slopes are given in the appendices. 



2 KLB picture for unbounded 2D fluids 

The evolution of the ensemble- averaged energy spectrum E(k), which repre- 
sents the energy density associated with the wavenumber k, is governed by 
(see Frisch [11] and Kraichnan [13]) 



—E(k) = T(k) - 2uk 2 E(k) + F(k). (1) 

(JjL 

Here T(k) and F(k) are, respectively, the ensemble-averaged energy transfer 
and energy input rate and v is the kinematic viscosity coefficient. Since waves 
of the same scale do not nonlinearly interact, T{k) is linear in the modal 
component corresponding to wavenumber k. Moreover, T{k) satisfies, by virtue 
of energy and enstrophy conservation, 

oo oo 

jT(k)dk = J k 2 T{k)dk = 0. (2) 
o o 

For a fluid in a doubly periodic domain, which for convenience we call a 
bounded fluid, the integral in (2) and elsewhere in this section should be 
replaced by a discrete sum over all wavenumbers. This constraint on T(k) 
imposes certain restrictions on its distribution and is thought to give rise to 
the dual cascade, believed to be a distinct feature of 2D turbulence. 

We multiply (1) by k 2 and integrate both the original and resulting equations 
over all wavenumbers, noting from (2) that the nonlinear terms drop out, to 
obtain evolution equations for the total energy density E = J °° E(k) dk and 
enstrophy density Z = J* °° k 2 E(k) dk, 
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d_ 

dt 
d_ 

dt 



IvP + r], 



2uZ + e, 



(3) 



(4) 



in terms of the energy and enstrophy injection rates e = / °° F(k) dk and rj = 
J °° k 2 F{k) dk, respectively, and the palinstrophy density P = / °° k A E{k) dk. 
We assume that the forcing is spectrally localized to a wavenumber interval 
[fcmin, fcmax] in the sense that 



This hypothesis is employed in [39]. This is a classical (although not exclusive) 
scenario for the KLB theory (cf. Kraichnan [13], p. 1421b; Pouquet et al. 
[28], p. 314; and Lesieur [16], p. 291), and is furthermore a common setup in 
numerical simulations of forced 2D turbulence (cf. Lilly [18]; Basdevant et al. 
[1]; Shepherd [33]). This assumption seems plausible for time-independent or 
white-noise forcing over [k min , k ma _ x ], given the assemble-averaged nature of e 
and rj. In particular, a monoscale forcing at a wavenumber s satisfies (5) with 
kmm = fcmax = s for each individual realization. Another example of such a 
forcing over [k m i n , k max ] is described in [39]; in each realization it yields time- 
independent energy and enstrophy injection rates e and rj such that rj = s 2 e, 
where s 2 is the mean of k 2 over [k m [ n ,k max \. A similar forcing was used by 
Shepherd [33] in a study of 2D turbulence in a large-scale zonal jet on the so- 
called beta-plane. In these examples, the characteristic forcing wavenumber s 
is constant in time. 

The dual cascade scenario can be best appreciated if one examines the evolu- 
tion equation 



obtained from (3) and (4). Here s is the forcing wavenumber yi]/e (for e > 
0), which according to hypothesis (5), must lie in the interval [k min , k inax }. 
If e = 0, then rj = and we take s to be any wavenumber in [k m i n , fc max ]. A 
direct enstrophy cascade requires that the characteristic enstrophy dissipation 



wavenumber y P/Z be much larger than the forcing wavenumber s; hence, the 
right-hand side of (6) must be positive. The positiveness of P — s 2 Z implies the 
positiveness of the left-hand side of (6) as well. If, in accord with the quasi- 
steady KLB theory, the total enstrophy reaches a steady state, we deduce 
that dE/dt must be positive. Equation (6) reflects the fact that the total 
energy must increase without limit, due to the inverse energy cascade toward 
wavenumber zero. (In a bounded domain, Tran and Shepherd [39] showed that 
P = s 2 Z and concluded from this that no direct cascade is possible.) 



< k 2 nhl t <rj< k 2 max e. 



(5) 




(6) 
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Equation (6) satisfies P/Z 3> s 2 not only for the KLB scaling A; -5 / 3 and A; -3 
but also for a rich variety of spectra. This condition only requires an energy 
spectrum shallower than k~ 5 for a sufficiently wide range of wavenumbers 
k > s, provided that the energy spectrum for k <C s is shallower than A; -3 . (In 
Section 4, it is shown that the quantity P — s 2 Z can be positive even if the 
energy spectrum for k > s is steeper than k~ 5 . This allows for the possibil- 
ity of an inverse energy cascade in the absence of a direct enstrophy cascade.) 
However, the KLB theory insists on the specific scalings A; -5 / 3 and A;~ 3 (with a 
logarithmic correction proposed by Kraichnan [13,14] and further investigated 
by Bowman [6]), respectively, for the energy and enstrophy inertial ranges. 
The /c~ 5 / 3 scaling is analogous to the Kolmogorov spectrum for 3D turbu- 
lence; the k~ 3 scaling implies that successive octaves in the enstrophy range 
contain equal amounts of enstrophy, so that the enstrophy grows logarithmi- 
cally with dissipation wavenumber k u . There exist a number of predictions for 
the numerical value of the enstrophy inertial range slope —(3 in the literature. 
Saffman [31] proposes {3 = 4, while Moffat [24] favours a slightly smaller value: 
(3 = 11/3. Sulem and Frisch [37] instead propose the upper bound (3 < 11/3. 

In the long-time limit, the KLB inverse energy cascade (or any inverse energy 
cascade with a spectral scaling shallower than A; -3 near k = 0) carries no 
enstrophy with it. Therefore, the enstrophy necessarily approaches an absolute 
equilibrium in a quasi-steady state. As a consequence, (6) reduces to 



s 2 ^E = 2u(P-s 2 Z). (7) 

This equation indicates a simple and interesting fact about the KLB theory 
for unbounded 2D turbulence: in a quasi-steady state (for which dZ/dt = 0), 
the strength of an inverse energy cascade (the rate of the energy growth), if 
realizable, is primarily determined by the rate 2vP of enstrophy dissipation. 



Besides the simultaneous conservation of energy and enstrophy, other essential 
features of 2D turbulence that underly the KLB theory are the scale-selectivity 
of the molecular viscosity and the unboundedness of the domain (in both 
directions). Together, they give rise to an infinite reservoir of energy in the 
vicinity of k = that allows for the possibility of the KLB inverse energy 
cascade (which contains no enstrophy if the spectrum near k = is shallower 
than A; -3 ). The theory breaks down when either the scale-selectivity of the 
dissipation or the unboundedness of the domain is absent. We demonstrate 
the former case and the semi-bounded case below. The case of a fluid in a 
doubly periodic domain is studied in [39] and will be reviewed in the next 
section. 

Consider (1) with the viscous dissipation term uk 2 replaced by a constant 
a > 0; this scale-neutral frictional dissipation is often called Ekman drag in 
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geophysical contexts. Equations (3) and (4) become 



^-E = —2aE + e, (8) 
J t Z = -2aZ + ,, (9) 

which, for bounded injection rates e and 77, implies that both the energy and 
enstrophy are bounded. This simple fact precludes the KLB type of inverse 
energy cascade. (As argued below, an inverse energy cascade dissipated by the 
friction at the large scales is not possible either.) Moreover, upon applying 
(5), we find 



j t (k 2 miQ E-Z)<-2a(kl, m E-Z), (10) 
j t (k 2 max E-Z)>-2a(k 2 max E-Z). (11) 
Hence, in the limit t — > 00 the following holds 



k 2 mhl E -Z<0< k 2 max E - Z, (12) 
or equivalently, 



&min < ^ < ^max- (13) 

Equation (13) implies that the redistribution of energy and enstrophy obeys 
exactly the same constraint as that imposed on the energy and enstrophy injec- 
tion rates. Now, the boundedness of energy (enstrophy) prohibits an infinitely 
wide range of wavenumbers k < k min (k > k max ) in which the energy spectrum 
can scale as k^ 1 {k~ 3 ), as this would imply a logarithmic divergence of the 
energy (enstrophy) as k — > (k — > oo). In fact, an energy spectrum k~ x [k" 3 ) 
or steeper (shallower) on the large (small) scales is inconsistent with (13). 
Spectra consistent with (13) require that the energy and enstrophy be pri- 
marily dissipated near the region of forcing; hence no inverse (direct) energy 
(enstrophy) cascade is possible. 

Besides the boundedness of the energy density in the present case, the spectral 
distribution of energy and enstrophy obeying (13) is profoundly different from 
that of the classical picture. The energy range is seen to be much shallower 
than A;~ 5 / 3 for a friction that acts uniformly on all scales. What seems curious 
is that the slopes of the enstrophy range in both cases do not differ by much. In 
fact, one could argue that a k~ 3 enstrophy-range spectrum with a logarithmic 
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correction is consistent with the constraint (13), provided that the spectrum in 
the energy range scales as k~ l (with a similar correction). Nevertheless, we see 
that the scale-selectivity of the molecular viscosity plays an important role in 
the dual cascade picture. Perhaps, instead of studying the 2D Navier-Stokes 
dual cascade, one could replace the molecular dissipation vk 2 in (1) by another 
scale-selective dissipation v^k 2 ^ (with /x > 0) and examine the realizability of 
a dual cascade in this hypothetical unbounded system. 

Finally, it is interesting to note that the KLB picture is not realizable for 
a 2D Navier-Stokes fluid confined to an unbounded channel with a periodic 
boundary condition in the across-channel direction (0 < y < L) and vanishing 
velocity (and derivatives thereof) at x — ±oo. This system is also furnished 
with the zero-mean flow condition 



L 

J u(x,y,t)dy = 0, (14) 
o 

where u(x,y,t) is the fluid velocity. To see why the dynamics of this system 
are incompatible with the KLB picture, consider the Poincare inequality for 
this domain: there exists a constant A > such that 



oo oo 

A J k»E(k)dk< J k 2+fl E(k)dk, (15) 
o o 

where \x > 0. The optimal value of the constant A depends on L. Now, the 
energy equation (3) and enstrophy equation (4) can be bounded via (15) as 
follows 



^-E < —2v\E + e, (16) 

^-Z<-2u\Z + rj. (17) 

These inequalities render the boundedness of both energy and enstrophy. The 
boundedness of the energy rules out the persistent inverse cascade needed for 
the positiveness of P — s 2 Z; rather, an absolute equilibrium is more plausible 
for this system. This result is quite physically reasonable: if one visualizes the 
inverse cascade as a result of the coalescence of same-sign vortices to form 
ever-larger ones, the system will tend toward equilibrium as the radii of the 
vortices approach the channel width L. This dynamical behaviour is observed 
by Rutgers [30] in an experiment of turbulence in a long channel, where, as 
equilibrium is approached, the vortices grow until the channel width is reached. 

Now, if an equilibrium is achieved, (6) necessarily reduces to 
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P-s 2 Z = 0. 



(18) 



This equation, which also applies to the case of bounded fluids in equilibrium 
considered in the next section, implies that the dissipation of enstrophy mainly 
occurs in the forcing region and admits much steeper spectra than the KLB 
spectrum. 

Remark 1. Some theoretical studies of the 2D Navier-Stokes equations in 
unbounded domains make an a priori assumption that the Poincare inequality 
holds (e.g. Rosa [29]; see also Temam [38], p. 307). Such studies automatically 
exclude the possibility of the KLB dynamics. 



3 Dynamics of bounded fluids in equilibrium 

In bounded systems, the energy is also in absolute equilibrium for arbitrary 
(but positive) viscosity coefficient u, so that (18) holds in equilibrium. This 
implies that P/Z = s 2 , ruling out the existence of an enstrophy inertial range, 
as argued in [39]. Moreover, the energy spectrum for k > k max is steeper 
than k~ 5 . This constraint clearly indicates a dramatic departure from the 
KLB theory. It is consistent with numerous numerical results, in which large- 
scale vortices, known as coherent structures, are observed (see for example 
Borue [5]; McWilliams [21,22]; Santangelo, Benzi, and Legras [32]; Smith and 
Yakhot [34,35]); these are often blamed for causing spectra steeper than those 
predicted by KLB. Although the mechanism behind these structures is not 
fully understood, we argue that it is the steepness of the spectrum (steeper 
than A; -5 ) that allows coherent structures to form, rather than the other way 
around. There is no need to invoke coherent structures to explain steep spec- 
tra: the steepness arises merely as a consequence of global conservation laws, 
molecular viscosity, and a spectrally localized forcing. As a matter of fact, 
a small-scale spectrum steeper than k~ 5 implies that the large scales carry 
virtually all of the system's enstrophy. Hence, dynamics exhibiting strong 
large-scale structures on a much weaker turbulent background of noise are 
consistent with [39]. It should be noted that the small-scale spectrum only 
needs to be steeper than k~ 3 for most of the enstrophy to reside in the large 
scales. Hence, large-scale structures may also be observed in simulations where 
a slope between —3 and —5 in the enstrophy range can be achieved, using a 
large-scale dissipation (see next section). Also, if the spectrum on the large 
scales is shallower than A; -3 , one may expect vortices comparable in size to 
the forcing scale to form, because most of the system's enstrophy is then dis- 
tributed in that spectral region. This has previously been noted by Paret and 
Tabeling [27]. 

In Section 4, we show that (18), which may be rewritten as 
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J> 2 - k 2 )k 2 E{k) = J2(k 2 - s 2 )k 2 E(k), 

k<s k>s 



(19) 



implies the spectral exponents satisfy (3 > 5 and a + [3 > 8, where —a {—0) 
is the asymptotic slope of the range of wavenumbers lower (higher) than the 
forcing wavenumber (in the usual case where the forcing wavenumber is no 
greater than the geometric mean of the integral and dissipation wavenumbers). 
For example, if the small-scale spectrum is approximately k~ 5 , then the large- 
scale spectrum should scale as k~ 3 . This would be consistent with the observed 
k~ 3 spectrum for the large-scale dynamics of the atmosphere (Lilly [19]). Now, 
a large-scale k~ 3 spectrum means that the enstrophy scales as /c -1 ; each octave 
in this range contains approximately the same amount of enstrophy. Therefore, 
the dissipation of energy is uniformly distributed among successive octaves in 
the energy range, so that no inverse energy cascade is possible. Nevertheless, a 
spectrum steeper than k~ 3 on the large scales is allowed by (19). This is more 
likely to occur if the small-scale spectrum is only marginally steeper than k~ 5 
(see Section 4). Unlike the KLB inverse cascade, which carries virtually all 
of the injected energy to ever-larger scales, this inverse cascade, if realizable, 
would only carry a fraction of the energy input to the largest scales. 

The dynamics of a bounded fluid in equilibrium is characteristically different 
from the quasi-steady KLB picture. There is an infinite energy reservoir at 
k = in unbounded systems that is forever available to collect the energy 
transfer; this feature is absent in bounded fluids (and in unbounded fluids 
satisfying the Poincare inequality). Hence, the departure should not come as 
a surprise. What seems ironic is that the enstrophy is only weakly dissipated 
in bounded fluids at high Reynolds numbers but strongly dissipated in the 
unbounded KLB case: in the bounded case the result P/Z = s 2 [39] implies 
that the enstrophy dissipation rate 2vPjZ becomes 2us 2 , while in the KLB 
theory it is approximately vk 2 / \n{k v / s), which is much greater than Ivs 2 since 
k v > s. 

There appears to be no simple generalization from the dynamics of a bounded 
fluid in absolute equilibrium to that of its unbounded classical counterpart and 
vice versa. The familiar reconciliation found in the literature is that the k~ 5 ^ 3 
range is modified or disrupted at the large scales when the inverse cascade 
reaches the lowest available wavenumber. If the KLB picture applies to non- 
equilibrium dynamics in a bounded system, v.z. before the inverse cascade 
gets reflected by the spectral boundary,^ then a dramatic adjustment of the 

2 In a finite system, an inverse energy cascade carries a non-negligible amount of 
enstrophy. If an inverse energy cascade carrying virtually all of the energy input 
rate e, reaches a wavenumber fc*, it would transfer enstrophy at the rate k 2 e to 
wavenumber k*. Hence, the ratio of the enstrophy accompanying the inverse energy 
cascade to the enstrophy input is approximately given by k 2 /s 2 , independent of the 
viscosity v. 
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spectrum has to occur as an equilibrium is approached (see Smith and Yakhot 
[35] for a discussion of the so-called finite-size effects). We have not yet inves- 
tigated in detail how this adjustment takes place; however, it seems plausible 
that as the inverse cascade hits the spectral boundary and gradually loses its 
strength (dE/dt — > 0), a substantial amount of the energy gets bounced back 
to the forcing scale. As there is little dissipation at the large scales, growth of 
the energy spectrum in the energy range is inevitable. This growth proceeds 
until the large-scale energy spectrum is sufficiently excited so that the enstro- 
phy dissipation occurs mainly in the vicinity of the forcing scale, whereupon 
a forced-dissipative equilibrium is reached. As dE/dt — >• 0, the enstrophy cas- 
cade (if initially present) ceases since the quantity P — s 2 Z on the right-hand 
side of (6) decreases to zero. As a consequence, a gradual steepening of the 
enstrophy-range spectrum takes place (whatever the spectral slope of the en- 
strophy range during the transient phase). In Section 4 and Appendix A, we 
establish that the sum of the steady-state spectral exponents in the energy 
and enstrophy ranges must asymptotically approach —8. 

Remark 2. It should be emphasized that the possibility of a direct enstrophy 
cascade in a bounded fluid during the non-equilibrium phase cannot be ruled 
out (not to say, however, that it actually occurs). For the case of a monoscale 
time-independent forcing, it is shown in Tran and Shepherd [39] that a direct 
enstrophy cascade is not realizable on average, whether the average be taken 
on a chaotic trajectory, limit cycle, or on the entire global attractor. This result 
leaves only the possibility of a direct enstrophy cascade in a neighborhood of 
the resulting monoscale stationary solution, on its unstable manifold. In this 
region one simultaneously has P — s 2 Z > and s 2 E — Z > (see Tran and 
Shepherd [39] and also Tran, Shepherd, and Cho [40]); the former inequality 
prevents one from ruling out a direct enstrophy cascade. More quantitative 
determination of the quantity P — s 2 Z (which is viscosity-dependent, see [40]) 
in this region will help assess the existence of the enstrophy cascade (and how 
this would depend on the viscosity) as a monoscale basic flow loses its stability. 

It is interesting to note that a k~ 3 spectrum (or slightly steeper), subject to 
experimental error, is observed in the laboratory experiments of Paret, Jullien, 
and Tabeling [25] , and Rutgers [30] . In these experiments, mechanical friction 
at the bottom and top boundaries (in particular with the noted in [30]) 

of the fluid could be sufficiently strong to outplay viscosity, so that Ekman drag 
alone is essentially responsible for dissipation. This might explain the observed 
spectra, according to the discussion in the previous section. Moreover, the 
analysis of [39] suggests that a combination of strong Ekman drag and weak 
viscosity allows for such a shallow spectrum to be realizable. Thus, there is no 
contradiction between the observed spectra and the predicted k~ 5 spectrum 
for a Navier-Stokes fluid in equilibrium. 
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4 Constraints on constant spectral slopes 



We now derive constraints on the spectral slopes for bounded fluids in equilib- 
rium and for a dual cascade in unbounded fluids, particularly for a persistent 
inverse energy cascade. We assume that the quasi-steady (steady) spectrum 
for the unbounded (bounded) case can be approximated by 



where a, b, a, j3 are constants, ko is the lowest wavenumber in the energy 
range (ko — > as t — > oo for the unbounded case), and k v is the highest 
wavenumber in the enstrophy range, beyond which the spectrum is supposed 
to be steeper than k' 13 . In Appendix A, we show that the arguments below can 
even be extended to the more realistic case, where the inert ial-range slopes 
depend on wavenumber. 

The quantity P — s 2 Z can then be estimated as 



J (k 2 - s 2 )k 2 E(k) dk>aj(k 2 - s 2 )k 2 ~ a dk + bj(k 2 - s 2 )k 2 ^ dk 

k k s 

1 1 

= aS 5 - a J (K 2 - l)K 2 - a dK + bs 5 -P J (\-K 2 )K^dK 

s/k v 

1 \ 

- j (1- K 2 )K 2 ~ a dK+ j (1 - K 2 )^' 6 dn , 
ko/s s/k v J 

(21) 



= as 




where the inequality results from dropping the spectral contribution beyond k v 
(which is considerable if f3 < 5), the second line is obtained by the respective 
changes of variables k = k/s and k — s/k in the two integrals on the right- 
hand-side of the first line, and the third line is obtained using the continuity 
relation as~ a = bs" 13 . 

In the bounded case the left-hand side of (21) vanishes. Therefore, 



j (l-K 2 )K 2 ~ a dK> J (1 - K 2 )^' 6 dK. (22) 
ko/s s/kv 

For a strong forcing at a relatively low wavenumber s, it is reasonable to 
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assume that k /s > s/k u . This requires 2 — a < {3 — 6, as the integrals in (22) 
decrease if the corresponding powers of k {(3 — 6 and 2 — a) increase. That is, 
a + {3 > 8. 



On the other hand, the convergence of the right-hand integral in (22), as 
s/k u — > 0, requires also that (3 > 5. This result was derived in [39], on the 
basis that the dissipation of enstrophy mainly occurs in the vicinity of the 
forcing scale. Now if f3 — 5 + 5, where it may be plausible that < 5 <C 1 for 
high Reynolds numbers, then a > 3 — 5. It thus seems possible to obtain a ~ 3. 
In the limits k /s — >■ and sjk v — > 0, as is usual for high- Reynolds-number 
turbulence, the inequality a + (3 > 8 approaches an equality. 



Remark 3. We caution that replacing the molecular viscosity vk 2 by a general 
viscosity v^k 2 ^ (fi > 0) in the previous argument leads to the result a + (3 > 
4 + 4/i. The significant dependence on \i of this constraint suggests that the 
introduction of a hyperviscosity could seriously alter the expected steady-state 
spectral slopes. 



We now derive a condition for a persistent inverse energy cascade in the un- 
bounded case. We assume a < 3, in accord with the realization of an inverse 
cascade toward wavenumber k — that carries no enstrophy with it. In the 
limits k /s — > and s/k u — > 0, the condition P — s 2 Z > is guaranteed if 
2 — a > (3 — 6 and a < 3, or equivalently, a + (3 < 8 and a < 3. These con- 
straints admit a variety of spectra for a quasi-steady state in which an inverse 
energy cascade to wavenumber zero, carrying with it virtually no enstrophy, 
and a direct enstrophy cascade to the dissipation wavenumber k u , carrying 
with it virtually no energy, are allowed. Note that the inverse cascade scenario 
cannot be ruled out even if (3 > 5, i.e, in the absence of a direct enstrophy 
cascade, provided that the inequality a + (3 < 8 holds. For the KLB energy- 
range spectrum k~ 5 ^ 3 , it is interesting to note that this condition requires only 
(3 < 19/3. We therefore suggest that an inverse energy cascade in the absence 
of a direct enstrophy cascade can be realizable for a wide range of (modest) 
Reynolds numbers. 



Remark 4. If the molecular viscosity vk 2 is replaced by a hyperviscosity 
v^k 2 ^ 1 for fi > 1 in an unbounded fluid, the condition for a persistent inverse 
energy cascade P — s 2 Z > is replaced by J^°{k 2 — s 2 )k 2 ^E(k) dk > 0. This 
leads to a + (3 < 4 + A/i. Of course, the condition a < 3 is required for a 
zero-enstrophy-carrying inverse energy cascade. 
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5 Large-scale dissipation 



In this section, we examine how a large-scale dissipation could be used to ob- 
tain a dual cascade, and in particular, a direct enstrophy cascade. Consider (1) 
in the bounded case, with a general dissipation: 



Y E ( k ) = T ( k ) ~ D(k)E(k) + F(k), (23) 

where D(k) is a non-negative function of k. Systems for which D(k) vanishes 
in the intermediate wavenumber range, including the forcing region, and for 
which the boundedness of energy is not guaranteed (Eyink [9]) have previously 
been studied in the literature. 

For a general D(k), (6) becomes 



s 2 jE-^Z^{k 2 -s 2 )D{k)E{k), (24) 
which, in equilibrium, reduces to the balance equation 



Y,(k 2 -s 2 )D(k)E(k) = 0. (25) 

k 

A slightly different form of this equation is derived in [39] . Equation (25) im- 
plies that the energy-range spectrum is related to the enstrophy-range spec- 
trum in an intimate manner. For a given D(k), an increase of the energy in one 
range requires an increase of the energy in the other. Thus, a steeper (shal- 
lower) energy-range spectrum corresponds to a shallower (steeper) enstrophy- 
range spectrum, for fixed E(s). Another obvious consequence of (18) is that a 
nontrivial equilibrium is not possible if D(k) vanishes for all k < s. 

We are interested in an expression for D(k) that retains the usual molecular 
viscosity and includes a large-scale dissipation. Thus, we consider D(k) = 
D e (k) + 2uk 2 , where D e (k) is a non-negative function of k. This dissipation 
includes the physically relevant case in which D^(k) is a positive constant 
representing friction from the planetary boundary layer in the geophysical 
context. Equation (25) then becomes 



2u(P - s 2 Z) = ]T(s 2 - k 2 )D e (k)E(k). (26) 

k 

The left-hand side of (26) is the familiar term due to viscosity, which would 
vanish in the absence of the large-scale dissipation. It can now become positive 
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since the right-hand side can be made positive in a variety of ways. Two 
popular forms of D^k) used in numerical simulations are the inverse viscosity 
Di(k) = Iv^k 2 ^ , for /i < 0, and mechanical friction restricted to the largest 
scales (say k < k e ): D^{k) oc H(k e — k), where H(ke — k) is the Heaviside step 
function (see Maltrud and Vallis [20], for example). To facilitate the formation 
of an enstrophy cascade, one might try to maximize the right-hand side of (26). 
However, since we have no a priori control over E{k) for different D^k), it is 
not known how to maximize the product Dp(k)E(k) for k < s. Nevertheless, 
restricting D e (k) to k < s by setting D £ (k) = for k > s is reasonable since 
any non-zero contribution to the right-hand side of (26) beyond s is negative. 

If a positive value for the right-hand side of (26) can be achieved, the quantity 
P—s 2 Z will be large for sufficiently small v. This may help break the constraint 
that the enstrophy range slope must be steeper than k~ 5 and allow for a direct 
enstrophy cascade. An inverse cascade should be realizable, as we expect most 
of the energy dissipation to occur at the large scales. Thus, a dual cascade is 
possible. However, this cannot be achieved without a cost, as De(k), which 
dissipates energy at the rate J2k De(k)E(k), also dissipates enstrophy on the 
large scales at the rate J2k D e (k)k 2 E(k). But 



k 2 Q £ D e {k)E(k) < £ D e (k)k 2 E(k), (27) 

k k 

where k is the lowest wavenumber corresponding to the system size. Hence, 
if J2k De{k)E(k) is comparable to the energy injection rate (which is ideally 
sought after in the spirit of the KLB theory), then the ratio of the enstrophy 
dissipation rate at the large scales to the enstrophy injection rate is greater 
than Amax- This fraction of the enstrophy dissipation at the large scales may 
be small, but not negligible. By allowing the enstrophy to be transferred to 
the large scales, a non-negligible amount of enstrophy may be trapped in the 
forcing region. If this is the case, then the spectrum in the forcing region has 
to adjust dramatically (since v is small and the large-scale dissipation is as- 
sumed to be weak around the forcing scale) to balance the trapped enstrophy. 
In Appendix B, we emphasize the difficulty of obtaining an enstrophy-range 
spectrum shallower than k~ 5 with a large-scale dissipation that is well sepa- 
rated from the forcing region. 

Numerical simulations of 2D turbulence can resolve up to a certain wavenum- 
ber, say kr- Therefore, there is always a finite amount of energy dissipation at 
the small scales. This is analogous to the dissipation of enstrophy at the large 
scales previously considered, due to the finite size of the domain. If Dh{k) 
represents the small-scale dissipation coefficient, the enstrophy and energy 
dissipations on the small scales are respectively given by J2k Dh{k)k 2 E{k) and 
J2kDh(k)E(k). These quantities satisfy 
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J2 D h {k)k 2 E{k) <k 2 T Y. D h {k)E{k). 



(28) 



k k 



Hence, if the dissipation of enstrophy by Dh(k) is comparable to the enstrophy 
injection rate (which is ideally sought after in the spirit of the KLB theory), 
then the ratio of the energy dissipation rate at the small scales to the energy 
injection rate is greater than k\^Jk\. This is true for any D h (k) > 0, including 
a hyperviscosity of arbitrary degree. 

Intuitively, if a large-scale dissipation is extended to s, one would expect it to 
absorb the reflected energy and keep the spectrum in the forcing region from 
growing as v — > 0. For a strong forcing and strong large-scale dissipation D e (k) 
(confined to k < s), it may be hypothesized that the value of the right-hand 
side of (26) is unaffected as v — > 0, given all else fixed. If this is the case, the 
quantity P — s 2 Z grows as v~ x and the ratio P/Z is given by 



This makes P/Z ^ oo as i/ ^ 0, a favorable limit for a direct enstrophy 
cascade, with a spectrum shallower than k~ 5 . However, the resulting cascade 
would not have the physical significance of the KLB theory since the direct 
enstrophy cascade (regardless of the spectral slope) might only be marginal, 
with a significant fraction of the enstrophy dissipated in the energy range 
(k < s) due to the strong large-scale dissipation, contrary to the classical 
theory. 



In this paper we have analysed the classical dual cascade theory of 2D tur- 
bulence in unbounded fluids formulated by Kraichnan [13,14], Leith [15], and 
Batchelor [2]. The main feature of the theory — the dual cascade — is contrasted 
to the behaviour of 2D turbulence in a region that satisfies the Poincare in- 
equality, such as a doubly periodic domain. It is shown that the dual cascade 
picture, if realizable, would strictly be an unbounded-system phenomenon. 
This important point is not adequately stressed and has often led to confu- 
sion in the literature. The familiar qualitative argument that the A; -5 / 3 range 
is modified or disrupted at the large scales when the inverse energy cascade 
reaches the largest available scale in a bounded system (assuming the applica- 
bility of the dual-cascade dynamics to the transient phase) is inadequate. Two- 
dimensional turbulence either in a doubly periodic domain or in an unbounded 
channel with a periodic boundary condition on the across- channel dimension 
does not behave in the manner predicted by KLB. In particular, the spectral 



^ = s 2 + ^U* 2 -k 2 Wk)E(k). 

k 



(29) 



6 CONCLUSION 
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slopes in such systems are found to satisfy (3 > 5 and a + (3 > 8, where —a 
{—(3) is the slope of the range of wavenumbers lower (higher) than the forcing 
wavenumber. This result is well supported by numerical simulations, which 
consistently find enstrophy-range spectra steeper than the KLB prediction 
and dynamically dominant large-scale structures (McWilliams [21,22]; San- 
tangelo, Benzi, and Legras [32]). It may even explain the observed large-scale 
k~ 3 spectrum in the atmosphere (Lilly [19], Boer and Shepherd [3]). 

We have shown that a dual cascade in unbounded fluids is possible if (3 satisfies 
3 < (3 < 5; this includes both the classical (3 = 3 scaling as an extreme limit, 
and also the theories of Saffman [31], Moffat [24], and Sulem and Frisch [37], 
which propose f3 = 4, (3 = 11/3, and (3 < 11/3, respectively. Moreover, in 
the absence of a direct enstrophy cascade {(3 > 5), an inverse energy cascade, 
corresponding to flows with low Reynolds numbers (which should be relatively 
easy to simulate) cannot be ruled out. 

The fundamental difference between bounded and unbounded fluids is that 
there is an infinite energy reservoir in the unbounded case, which allows a per- 
sistent inverse energy cascade to ever-larger scales to form, so that the energy 
eventually evades viscous dissipation altogether. Provided that the spectrum 
near k = is shallower than k~ 3 , the inverse cascade asymptotically carries 
no enstrophy. This luxury is a consequence of both the unboundedness of the 
domain (in both directions) and the scale-selectivity of the molecular viscosity. 
In addition to the simultaneous conservation of energy and enstrophy, these 
properties constitute the basic building blocks of the KLB theory. Another 
important hypothesis is the existence of a quasi-steady state. As long as an 
inverse cascade is realizable and a quasi-steady state can be established in 
an unbounded system, the cascade dynamics are fundamentally distinct from 
what occurs in a bounded fluid in equilibrium. There appears to be no sound 
basis for extending the results from one case to the other. Of course, in an 
unbounded fluid it is quite possible for an inverse energy cascade to exist with- 
out a corresponding direct enstrophy cascade; in this case, there might then 
be certain similarities between the dynamics of the bounded and unbounded 
systems. 

The dissipation operator plays an important role in the spectral distribution of 
energy. This is especially apparent in the balance equation (25) for a fluid in a 
doubly periodic domain: the product of the energy spectrum and the spectral 
dissipation function in the energy and enstrophy ranges are intimately related. 
This information is important for numerical 2D turbulence simulations, where 
various dissipation mechanisms are employed: it should help researchers rule 
out certain spectral distributions and anticipate possible outcomes for a given 
dissipation mechanism. Finally, we showed that a large-scale dissipation could 
give rise to a direct enstrophy cascade since the quantity P — s 2 Z could grow 
as until P/Z > s 2 . 
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A Constraints on general spectral slopes 



Strictly speaking, the spectrum (20) is too simplistic; actual spectral slopes will 
tend to vary monotonically with wavenumber (particularly in the enstrophy 
range, as one approaches the onset of the dissipation range). The arguments 
of Section 4 can be readily extended to more general spectra. Normalizing all 
wavenumbers so that s = 1, we express the energy spectrum as 



' k~ a ^ if k < k < 1, 
E(k) = a{ k-PW if i < k < k u 

k~^ if fci < k < k T , 



(A.l) 



with /3(fci) = 7(fci), where k is the lower spectral cutoff wavenumber (de- 
termined by the domain size), k\ is the highest wavenumber in the enstro- 
phy range, and fcy is the highest retained (truncation) wavenumber. Equa- 
tion (18), or equivalently (25), then appears as 



1 fci kx 

J (I- k 2 )k 2 - a ^dk = J(k 2 ~ l)k 2 ~ m dk + J{k 2 - l)k 2 -^ k) dk. (A.2) 
fco i fci 

The change of variable k — 1/fc in the integrals on the right-hand side yields 



i i 

J(l- K 2 )K 2 ~ a ^dK = J (1 - K 2 )K^ K) ~ 6 dK + €„, (A.3) 
fco 1/fci 

where e v = fij^(l — /c 2 )/^ 1 /") -6 dk > 0. For bounded turbulence, we restrict 
our attention to the usual case where fci > 1/fco- Since Ji/ kl (l — K 2 )K e dK is 
a strictly decreasing function of 9, we find that the maximum slopes a = 
supa(fc) and (3 = sup/3(fc) satisfy 

k k 

i i 

J (\-K 2 )K 2 ~«dK> J (1 
1/fci 1/fci 

1/fcl 

Hence 2 — a < (3 — 6; that is, a + (3 > 8. 



l 

- K 2 )K 2 ~ a ^ dK> J (I- K 2 )K 2 ' a ^ dK 

fco 

1 

- k V (1/, °~ 6 ^ > / (! - K V~ 6 

i/fci 
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One can also obtain estimates for the minimum slopes a = inf a(k) and (3 = 

k — 

inf (3(k), assuming 'j(k) > (3: 

k — 

1 1 1 

|(l-«V~-^< j(l-K 2 )K 2 ~ a ^dK= j (l-K 2 )K^- 6 dK + e u 
ko ko 1/ki 

1 

< J (l-K 2 )t£- % dK. (A.4) 

l/k T 

In the asymptotic limit as ko — > and kr — * oo, we deduce a + (3 < 8. 

It is instructive to specialize these results to the case of constant slopes, where 
a = a = a and (3 — (3 — (3 (e.g. if k x <C k u ). For a bounded fluid satisfying 
k\ > 1/ko, we then see that lim^^o ot + (3 = 8. 



B Spectral slopes for a system with an inverse viscosity 

It was suggested in Section 4 that a large-scale dissipation, well separated from 
the forcing scale, may not give rise to the desired direct enstrophy cascade and 
the corresponding k~ 3 spectrum. To demonstrate this point, we consider the 
special form D(k) = 2v'k~ 2 + 2uk 2 , with v's~ 2 = us 2 , so that at the forcing 
scale, the inverse viscosity and the usual molecular viscosity have the same 
strength. Equation (25) becomes 

J2(k 2 -s 2 )(u'k~ 2 + uk 2 )E(k) = 0. (B.l) 

k 

Assume the spectral scaling (20) and follow the steps leading to (22), we obtain 

i i 
J (1 - k 2 ){k~ 2 + K 2 )K~ a dn> J (1 - k 2 )(k 2 + k- 2 )k?- 4 dn. (B.2) 

ko/s s/k v 

In the case ko/s > sjk v it follows that a+(3 > 4. In the limit s/k u — > 0, we still 
require f3 > 5. For the KLB enstrophy-range spectrum to be realizable it is 
necessary that a > 1; moreover, s/k u cannot be much smaller than ko/s. The 
former condition does not seem to be plausible in the presence of an inverse 
viscosity, while the latter condition requires an unphysically narrow enstrophy 
range for a forcing at relatively small wavenumbers. 
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